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ABSTRACT 

A theoretical study of the temperature dependence of the rate con- 
stants for model reacting systems has been carried out in an attempt to 
understand some recent experimental measurements which imply the existence 
of negative activation energies. A collision theory model and classical 
trajectory calculations are \ised to demonstrate that the reaction proba- 
bility can vary Inversely with collision energy for bimolecular reactions 
occurring on attractive potential energy surfaces. However, this is not 
a sufficient condition to ensure that the rate constant has a negative 
temperature dependence. On the basis of these calculations it seems un- 
likely that a true bimolecular reaction between neutral molecules will 
have a negative activation energy. 
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INTRODUCTION 


In recent years experimental rate constant measurements of some 
apparent blmolecular reactions have resulted in the assignment of negative 
activation energies. These so-called negative activation energies 
arise from the fitting of the experimental data to the Arrhenius 
equation: 


K(T) - A exp(-E^/kT). (1) 

The pre-exponential factor A and activation energy are constants deter- 

mined in the fit. Several hypotheses can be advanced to explain the 
measurements: (I) The reactions proceed along attractive potential 

energy surfaces without energy barriers on which slower collisions are 
more reactive than faster ones;^ (II) The reactions proceed via the 
formation of metastable Intermediates and, thus, are not true blmolecular 
reactions; and (III) The experimental data contain systematic errors. 

The present paper is an examination of the first hypothesis. This 
examination is made in terms of collision theory and classical trajectory 
calculations using model potential energy surfaces. We first examine 
whether slower collisions can be more reactive than faster ones, and 
second, whether this effect is sufficient to result in negative activa- 
tion energies. 

Rate data for some reactions for which negative activation energies 
have been reported are given in Table I. Note that most of the data have 
been reported with extremely small uncertainty limits. It is only in the 
last few years that such precise measurements of very fast gas phase 
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reaction rates have been possible. This advance In chemical kinetics 
research has arisen largely because of the concern over the detrimental 
effect of anthropogenic halocarbons on the earth's ozone layer. Comparl- 
sons of the compendia of rate constants for stratospheric reactions from 
1974 (Ref, 7) and 1977 (Ref. 8) reveal the extent and improved precision 
of the new body of chemical kinetics data. Of particular Interest are 
the measurements of activation energies over the temperature range 
200-500 K for some very fast reactions. Previously, uncertainties of 
0.5 to 1.0 kcal/mole existed In most measured activation energies and 
the small negative activation energies of interest In the present work 
were lost In the experimental uncertainty. 

The reactions listed in Table I can be analyzed in terms of the 
three possible explanations given above. A collision complex mechanism 
can be ruled out for reaction (Rl) since the intermediate CtOO species 
is very weakly bound with respect to the products and could have only a 
very short lifetime. ^ Therefore, the result of Zahnlser and Kaufman,^ 
if preferred over the experimental data of Clyne and Nip^ or the calculated 
results of Jaffe,® would In^ily the existence of a negative activation 
energy in a true bimolecular reaction. On the other hand, (R2) probably 
does proceed via a complicated mechanism involving at least two metastable 
intermediates.^ Cox' has reported preliminary rate constants for reaction 
(R3). Since it is a hydrogen transfer reaction, a mechanism involving an 
intermediate complex seems unlikely. However, the NH 2 -WO reaction (R4) 
undoubtedly does involve a complicated mechanism and the overall rate 

•a 

expression cannot be taken as being that of a bimolecular process. 

Davis et al.** have reported negative activation energies for the reactions 
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of atonic oxygen with tetranethylethylene (TME) ai^ (els - 2 - butene), 
reactions R5 and R6, respectively. They suggested** that these data 
could be explained by a kinetic nodel In which the energy dependence of 
the reaction cross section Is proportional to a delta function centered 
on the threshold energy. 

Thus, the obser^d rate constants for (R2) and (R4) can be explained 
In terns of complex mechanisms (hypothesis 11). The rest of the data 
presented In Table 1, If correct, must be explained In terms of a 
collision theory nodel In which slower collisions are more likely to react 
than faster ones (hypothesis 1). In the present work, we examine such a 
model in detail to determine If it Is capable of reproducing the observed 
reaction rate data. 

In the next section, we define a suitable collision theory model and 
demonstrate the conditions that must be met for a negative activation 
energy to be realized. In Sec. Ill, we describe four model potential 
energy surfaces based on the CIO •¥ 0 reaction (Rl) . These surfaces are 
used for classical trajectory calculations of reaction cross sections and 
rate constants which are described In Sec, IV. The resulting cross 
sections and rate constants are used to determine the reliability of the 
collision theory model and to determine whether negative activation 
energies should be computed for physically reasonable systems. The con- 
clusions of this study are given In Sec. V. 


II. ACTIVATION ENERGIES IN COLLISION THEORY 

The Arrhenius activation energy of a chemical reaction Is actually 
an empirical concept that arises from the fitcltig of experimental rate 
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constant data. In the cose of a system with a single Internal state 
(v,J), the activation energy has been shown to be equal to the average 
translational energy for reactive collisions minus the average trans- 
lational energy for all collisions. For the general case in which many 
reactant internal states are populated, one can show that is the 

difference in average total energy between reactive collisions and all 
collisions. Thus, negative activation energies can occur if enhanced 
reagent energy inhibits rather than promotes reaction. In this section, 
we construct a collision theory model to describe blmolecular reaction 
kinetics and examine the resulting activation energy. 

The activation energy is formally defined by inverting the Arrhenius 
rate constant expression [Eq. (1)]: 


.dlnK(T) 

“ ciTTr 


( 2 ) 


This simple expression Is generally not flexible enough to represent the 
variation of the rate constant over a large temperature range and a three- 
parameter expression is often used instead 


K(T) - BT** exp(-Ej^/kT) , (3) 

where B, q, and are varied to fit the experimental data. Substituting 
Eq. (3) Into Eq. (2), we have 


E. - qkT + E? . 
A A 


(A) 


The first term in Eq. (4) arises from the t«nperature dependence of the 
preexponential factor of the rate constant and is thought to be responsible 
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for the observed negative activation energies of the reactions listed in 


Table I. 

In order to examine the physical significance of q and wt con- 
sider the collision theory rate constant for a thermal atom-diatom 
exchange reaction (A+BC-+-ABrK:) 



/ flO 

S^(E,v,J)exp(-E/kT)EdE, 


where E and Q , are the BC vibration-rotation energy and partition 
v,J ^v,J 

function, respectively, ^ is the atom-diatom reduced mass, and 
is the total reaction cross section for BC in the (v,J) internal state 
ivud collision energy E. The reaction cross section depends on the 
Interatomic forces and must be determined by quantum or classical mechanical 
scattering calculations. For the purpose of this discussion, however, we 
define a simplified model cross section as follows: 


S^(E) - 0 , E < E^ 

S^(E) - aE^, E > E^ , 


( 6 ) 


where the paramenters a and f are independent of temperature and 
internal state (v,J) and is the threshold energy. Substituting (6) 

into (5) and defining X - E/kT, we have 


K(T) 



(kT) 


f+l/2 




X^+^e-*dX 


(7) 


6 



The integral In Eq. (7) Is the incomplete gamma function r(£+2,X^).^^ 
The activation energy defined in Eq. (7) can be shown to be 


E 


A 


(f + l/2)kT + 



f+i 

e;tp(-E^/kT) . 


( 8 ) 


Thus, q « f + 1/2 and the second term in Eq. (8) represents a temperature 

dependent E^. The special case £ * 0 is equivalent to the hard sphere 

1 /2 

cross section. This results in q ■ 1/2 and represents the T 
"collision frequency" factor coimiK)nly used in rate constant expression. 

Since > 0, E^ cannot be negative in this model [see Eq. (8)] 
and a negative activation energy can only arise if £ < -1/2. For most 
chemical reactions, however, E^ > kT and the second term in Eq. (8) 
usually dominates. The threshold energy E^ is approximately equal to 
the potential energy barrier which separates reactants and products. ^ 
It is believed that there are no such barriers for the systems listed in 
Table I, however, in which case E^ * 0. Setting » 0, the model rate 
constant for these reactions is given by 


K(T) 



(kT)^'*’^^ ^ r(f + 1/2) 


(9) 


Equation (9) is finite for f > -2 and E^ ■ (£ + l/2)kT is negative for 
-2 < f < - 1/2. Thus, the cross section must have a strong Inverse 
dependence on energy before a negative activation energy can be achieved. 
Since E* is physically similar to an activation energy one can say that 
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the rate constant given by Eq. (9) has a "zero activation energy" but a 
nonzero "tmperature dependence" equal to i ■¥ If 2. 

In the above model, the case of ■ 0 and f < 0 corresponds to 
the situation we are examining in the present study. For this case, slower 
collisions have larger cross sections; the larger cross sections mean they 
have greater probability of reacting than faster ones. However, only if 
£ < -1/2 will a negative tenq>crature dependence or Arrhenius activation 
energy result. 

III. MODEL POTENTIAL ENERGY SURFACES 

The calculation of reaction cross sections and rate constants re- 
quires knowledge of the potential energy surface of the system ur Jer study. 
Ab initio calculations of these surfaces are quite difficult and unneces- 
sary for the purpose of the present work since we only want to study a 
general phenomenon In chemical kinetics. Thus, we have constructed four 
empirical potential energy surfaces for a generalized reaction A+BC AB+C. 
These surfaces have been modeled loosely after reaction (1) and have a 
potential well corresponding to the ClOO equilibrium geomet’^’^® which is 
63.6 kcal/mole lower in energy than the reactants asymptote and 7.8 kcal/ 
mole lower than the product asynq)tote. The potential energy is given by 
the sum of Morse potentials for the AB and BC bonds and a trlatomlc 
double minimum bending potential which Is switched off as A aiul BC 
separate. 
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This bending potential takes the form 




f(Rab> 




c + (e^-1.)' 


+D 


( 10 ) 


where F(R.„) Is the switching function 
Ad 

F(R^) - 1 . < P - 6p , 

6p (R»n“P) 

F(R._) - 1/2 ^ . p - 8p < R.„ < p + 5 p ; (11) 

“ «p)2 + (R^-p)2 “ 


F(R^) - 0 


, R^ > p + 6p . 


Both the bending potential and the switching function have been discussed 
previously.® The bending potential is parameterized to reproduce 

the CdOO bond angle (110") and bending force constant.^® The choice of 
P and iSp determine the range of R^ over which the bending potential 


is switched on 

during A ■»- BC 

collisions. 

If P 

is small then 

A + BC can 

approach close 

enough to feel 

substantial 

A - B 

attraction for 

any BC 


orientation. This case corresponds to the most attractive potential sur- 
face. If p is large the triatomic bending potential will dominate the 
kinetics and only collisions occurring over a small range of 9^^ can 
be reactive. The width of this range of bond angles will be greater for 
collisions with higher collision energy, thereby introducing a positive 
temperature dependence to the reacc^ .i rate constant. Four potential sur- 
faces differing only in the choice of p and 6p were used in the present 
study. In addition an exponential repulsion term for R^^, 
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5.0 exp[-4(R^^-1.6)] , 


( 12 ) 


REP 


where distances are in A and energies in kcal/aole, was added to the 
potential energy expression to pr eut R^^ fron going to zero during 
A + BC collisions. The parameters used in constructing the potential 
energy surface are given in Table II. 

The potential energy profile along a schematic reaction coordinate 
is shown in Pig. 1 along with the four ranges of R._ spanned by the 
switching functions considered in the present study. Figures 2 and 3 
ere cjCs through these potential energy surfaces 

extreme cases (I and IV) at R^^ - 1.83A. Note that in Case I (Fig. 2) 

A and BC can approach to R^ < 2.5A without feeling any appreciable 

bending potential while in Case IV (Fig. 3) the V_ extends to R._ - aA. 

o Afi 

The reactants are still at the asymptotic energy when R._ « 4A but 
they are more than 2 kcal/mole down into the trlatomic potential well when 
R^p = 2.5A. The rate constants for reactions occurring on these potential 
surfaces should exhibit strikingly different temperature dependences. 


IV. CLASSICAL TRAJECTORY CALCULATIONS 

Classical trajectory calculations have been performed t ■ determine 
total cross sections and rate constants for reactions occurring on the 
icodel potential energy surfaces. These calculations will assess the 
validity of our premise and the collision theory model presented in Sec. II. 

The quasiclasslcal trajectory method^ was used to compute the 
cross sections as a £u:ctlon of collision energy for the most attractive 
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model potential surface (Case I) at collision energies between 0.1 a»*'^ 
l.O kcal/fflole and BC Internal energy state (v,J) * (0,15). These 
conditions correspond to 0.624 kcal/mole vibrational energy and 0.308 
kcal-mole rotational energy in BC. At each collision energy the cross 
section was based on the calculation of 250 trajectories with a maximum 
impact parameter of 5X. The resulting data, which are shown in Fig. 4, 
are well fit by ■ 27.0/E^/** X^. Analysis of these calculations showed 
that slower collisions can be reactive over a larger range of impact 
parameters than faster ones. For example the largest observed Impact 
parameter for a reactive trajectory decreased from 5.0 X to 4.5 X as E 
increased 0.1 to 1.0 kcal/mole. This substantiates the first part of our 
premise: namely, that the long-range attraction is nmre effective in 
promoting reaction at lower collision energies. Both the effective ra.tge 
of impact paran«ters and the total reactive cross sections are greater 
for slower collisions than faster ones. t now remains to be seen whether 
this effect can lead to a race constant with a negative temperature 
dependence. 

The phase space trajectory method^ was used to calculate thermal 
rate constants over the temperature range 200 K to 350 K for reactions 
occurring on eac.i of the four potential energy surf.'ces. Each rate 
constant is baseu on a sample of 500 trajectories and has a standard error 
of 3.2 tc 6.2% owing to the statistical sampling method used. 

Arrhenius plots of these data are shown in Fig. 5. It can be seen 
that the rate constants for Cases I to III exhibit almost no temperature 
dependence while the data for Case IV have a definite negative slope. 

Least squares fits to these data for log K(l/T) yield activation energies 
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of 0.05, 0.00, and 0.01 kcal/mole for Cases 1 to III, respectively. 

Thus these rate constants exhibit a small positive or zero temperature 
dependence. For Case IV we obtain an activation energy of 0.28 kcal/mole 
by a similar procedure. If we assume the rate constant expression given 
by Eq. (3) with » 0 then K where q is determined by least 

squares fitting to be 0.09, 0.00, 0.03, and 0.54 for the four cases, 
respectively. The calculated cross sections for Case I also indicate 
that E^ ~ E^ * 0. Using these data in the collision theory model, we 
arrive at the prediction that q * 0.25. 

The observed temperature dependences would indicate that ~ E”°*^ 
for Cases I to III according to the collision theory model presented in 
Sec. II. These results indicate that changes in the range of the bending 

potential do not seem to affect the temperature dependence of the rate 

constant or the form of the reaction cross section as long as the range 
of Vg does not exceed the range of the A-B attraction. For Case IV 
we have ~ In this case the colliding reactants feel the bend- 
ing potential before the onset of A-B attraction. 

The agreement is quite satisfactory considering that the classical 
trajectory rate constants include contributions from a thermal distribution 
of reactant vibration-rotation levels while the cross sections were 
determined for a single arbitrary internal energy state. The results of 
the trajectory calculations, then, support the conclusions of the collision 
theory model. The probability of reaction can vary inversely with 
collision energy for systems with attractive potential energy surfaces. 
However, this is not a sufficient condition to ensure that the rate con- 
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stant will exhibit a negative tei^>erature dependence for the potential 
energy surfaces considered. 


V. CONCLUSIONS 

The restilts of tl« total reaction cross section and rate constant 
calculations performed in this study demonstrate that slower collisions 
can lead to enhanced reaction probabilities. However, this fact by itself 
is not sufficient to cause a negative teuperature dependence in the rate 
constant. Indeed, none of the rate constants conQ>uted in this study 
exhibited such a temperature dependence. In addition, the simple collision 
theory model presented here, which examines the relationships between cross 
sections and activation energy, seems to qualitatively account for all 
r.'ese results. 

Thus, we conclude that the cross section must have a strong inverse 
dependence on collision energy in order for the rate constant to show a 
negative temperature dependence. Undoubtedly, the strength of the long- 
range force between reactants is an important factor. In ion-molecule 
reactions, where the intermolecular forces are stronger, the cross section 
will probably have a stronger inverse dependence on collision energy and 
negative activation energies will result. However, for neutral -neutral 
systems it seem unlikely that a true blmolecular rate constant will have 
a negative temperature dependence. 
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TABLE I, RATE CONSTANT DATA FOR REACTIONS WITH NEGATIVE ACTIVATION ENERGIES. 
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^fined according to K(T) 


TABLE II. POTENTIAL ENERGY SURFACE PARAMETERS. 



D (kcal/fflole) 
e 

R (X) 
6 


n 

Morse parameters 




AB 

63.59 

1.23 

3.313 

BC 

7.76 

1.83 

3.458 

b 

Vg parameters 




A 

31.1327 kcal/mole 



B 

192.0009 kcal/mole 



C 

0.9908 



D 

-123.7830 kcal/mole 



Case I 

II 

III 

IV 

p ,X 2.0 

2.5 

2.5 

2.5 

6p,A 0.5 

0.5 

1.0 

1.5 


^Parameters for '"£00 from Arkell and Schwager.^^ 
^From Jaffa.® 
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FIGURE CAPTIONS 


FIG. 1. Schematic reaction path for the reaction A + BC AB + C. The 
limits of the switching regions (p ± 6p) of the attenuation factor 
used In the bendlug potential for the four model potential surfaces 
are shown. 

FIG. 2. Contours of equal potential energy In kcal/mole as a function 

of R._ and 6.„_ for Case I. The limits of the switching region, 

Ao AoL 

R p ± iSp, are shown by dashed lines. Regions of the potential 
surface corresponding to reactants (A ■¥ BC) and Intermediate complex 
(ABC) are shown. The BC bond length Is held fixed at Its equilib- 
rium value of 1.83 A. 

FIG. 3. Contours of equal potential energy In kcal/mole as a function 

of R^_ and for Case IV. The entire plot falls within the 

Ad AdL 

limits of the switching region which is centered at R^ ■ p ■ 2.5 A. 
The BC bond length Is held fixed at Its equilibrium value of 1.83 A. 
FIG. 4. Quasl-classical trajectory reaction cross sections versus collision 
energy for Case I with v ■ 0, J » 15, The error bars reflect the 
statistical sampling error In the trajectory calculation. The solid 
curve Is the least squares fit to the cross section data. ~ 
is shown for comparison. This Is the relation which would give no 
temperature dependence according to the collision theory model. 

FIG. 5. Arrhenius plot of the phase space trajectory rate constant data 

for the four cases considered. The error bars reflect the statistical 
sampling error In the trajectory calculations. The curves are least 
squares fits to the rate constant data. 
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FIG. 1. Schematic reaction path for the reaction A + BC -► AB + C. The 
limits of the switching regions (p ± 6p) of the attenuation factor 
used In the bending potential for the four model potential surfaces 
are shown. 
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V ^ ABC) for Case I 

Rbc” 



2^ ABC, deg 

FIG. 2. Contours of equal potential energy In kcal/mole as a function 
of and O^bc limits of the switching region, 

Rab " p ± 6p, are shown by dashed lines. Regions of the potential 
surface corresponding to reactants (A + BC) and Intermediate complex 
(ABC) are shown. The BC bond length Is held fixed at Its equlllb<' 
rlum value of 1.83 A. 
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V (Rad, ^ ABC) for Case IV 


A + BC CASE IV 


1/5 5 \ 0 0 / 5/10 



120 180 2m 

K ABC, deg 


300 


FIG. 3. C-^ntours of equal potently] energy in kcal/mole as a function 
of R/tB ®ABC Case IV. The entire plot falls within the 

limits of the switching region which Is centered at ■ p * 2.5 X. 

The BC bond length Is held flx^d at Its equilibrium value of 1.83 X. 





E, Kcal/mole 

FIG. 4. Quasl-classlcal trajectory reaction cross sections versus collision 
energy for Case I with v ■ 0, J ■ 15. The error bars reflect the 
statistical sampling error In the trajectory .alculatlon. The solid 
curve is the least squares fit to the cross section data. 
is shown for comp'^rison. This Is the relation which would give no 
temperature dependence according to the collision theory model. 
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